Exploring the sub-eV neutrino mass range with supernova neutrinos 
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A new method to study the effects of neutrino masses on a supernova neutrino signal is proposed. 
The method relies exclusively on the analysis of the full statistics of neutrino events, it is independent 
of astrophysical assumptions, and does not require the observation of any additional phenomenon 
to trace possible delays in the neutrino arrival times. The sensitivity of the method to the sub-eV 
neutrino mass range, defined as the capability of disentangling at 95% c.l. the case m„ = 1 eV 
from m„ = 0, is tested by analyzing a set of synthetic neutrino samples modeled according to the 
signal that could be detected at SuperKamiokande. For a supernova at the Galactic center success 
is achieved in more than 50% of the cases. It is argued that a future Galactic supernova yielding 
several thousands of inverse (3 decays might provide enough information to explore a neutrino mass 
range somewhat below 1 eV. 

PACS numbers: 14.60.Pq, 97.60.Bw 



During the past few years, the large amount of data 
collected by solar 0] and atmospheric Q neutrino ex- 
periments provided strong evidences for non vanish- 
ing neutrino masses. The recent KamLAND result 
|3j on the depletion of the neutrino flux from nuclear 
power plants in Japan, gave a final confirmation of 
this picture. Since in the Standard Model (SM) of 
particle physics all neutrino species are massless, this 
constitutes the first direct evidence for new physics, 
and provides important informations for developing 
theories beyond the SM. However, to date all the ev- 
idences for neutrino masses come from oscillation ex- 
periments, that are only sensitive to mass square dif- 
ferences and cannot give any informations on single 
mass values. The importance of measuring the abso- 
lute value of neutrino masses cannot be understated. 
It is presently being addressed by means of a remark- 
ably large number of different approaches, ranging 
from laboratory experiments to a plethora of meth- 
ods that relay on cosmological considerations. Re- 
cent reviews can be found in Q. From the study of 
the end-point of the electron spectrum in tritium (3- 
decay, laboratory experiments have been able to set 
the limit m„ e < 2.2 eV This is already close to 
the sensitivity limit of on-going experiments. If neu- 
trinos are Majorana particles, the non observation of 
neutrinoless double j3 decay can constrain a partic- 
ular combination of the three neutrino masses. In- 
terpretation of these experimental results is difficult 
due to large theoretical uncertainties related to nu- 
clear matrix elements calculations. This is reflected in 
a model dependent limit m c J f < 0.35 - 1.24 eV 0,0. 
A tight bound m Vi < 0.7 eV was recently set by 
the WMAP Collaboration by combining measure- 
ments of cosmic microwave background anisotropies 
with data from bright galaxies redshift surveys @ and 
other cosmological data. However, this limit becomes 
much looser if the set of assumptions on which it re- 



lies is relaxed (see |9j for discussions of this point). 
Therefore it is quite important to keep looking for al- 
ternative ways to measure the neutrino masses. Model 
independent limits can be more reliably established by 
combining complementary experimental informations 
and, needless to say, it would be of utmost impor- 
tance if the mass values could be measured by means 
of more than one independent method. 

Already long time ago it was realized that Super- 
nova (SN) neutrinos can provide valuable informations 
on the neutrino masses [l0|. The basic idea relies on 
the time-of-flight delay At that a neutrino of mass 
m v and energy E v traveling a distance L would suffer 
with respect to a massless particle: 
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where for ultra-relativistic neutrinos we have used 
l/v — E u /p v ~ l + ml/2E^. The dispersion in the ar- 
rival time of about twenty electron anti-neutrinos from 
supernova SN1987A was used in the pas t to set the 
model independent limit mp e < 30 eV |ll|. This limit 
can become significantly tighter under some specific 
assumptions. For example a recent detailed reanaly- 
sis obtained, within the SN delayed explosion scenario, 
the limit m Pe < 5.7 eV 12]. 

Since SN1987A, several efforts have been carried out 
to improve the sensitivity of the method, while await- 
ing for the next explosion within our Galaxy. Often, 
these approaches rely on "timing" events related to 
the collapse of the star core, that are used as bench- 
marks for measuring the neutrino delays. The emis- 
sion of gravitational waves |l3l Il4j , the v e neutroniza- 
tion burst 1141. the initial steep raise of the neutrino 
luminosity [l5j, and the abrupt interruption of the 
neutrino signal due to a further collapse into a black 
hole have been used to this aim. However, there 
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are some drawbacks to these methods: firstly only 
neutrinos with arrival time close to the benchmarks 
are used, and they represent only a small fraction of 
the total; secondly the observation of the benchmark 
events is not always certain, and in any case some 
model dependence on the details of the SN explosion 
is generally introduced. 

The method we want to propose is free from these 
drawbacks: it relies only on the measurement of the 
neutrinos energies and arrival times, and it uses the 
full statistics of the detected signal, thus allowing to 
extract the maximum of information. It is also re- 
markably independent of particular astrophysical as- 
sumptions, since no use is made of benchmarks events. 
The basic idea is the following: in the idealized case of 
vanishing experimental errors in the determination of 
the neutrino energies and arrival times, and assuming 
an arbitrarily large statistics and a perfectly black- 
body neutrino spectrum, one could use the events with 
energy above some suitable value E* (to suppress the 
mass effects) to reconstruct very precisely the evolu- 
tion in time of the neutrino flux and spectrum. Once 
the time dependence of the full signal is pinned down, 
the only parameter left to reconcile the time distribu- 
tion of the low energy neutrinos with the high energy 
part of the signal would be the neutrino mass, that 
could then be nailed to its true value. Of course, none 
of the previous conditions is actually fulfilled. In wa- 
ter Cerenkov detectors as SuperKamiokande (SK) the 
uncertainty on relative timings is negligible; however, 
the errors in the energy measurements are important 
and must be properly taken in to account. The statis- 
tics is large but finite, and this not only represents a 
source of statistical uncertainty, but also implies an 
upper limit on the useful values of E*. Finally, the 
SN neutrino spectrum is not perfectly thermal |l7j . 
Nevertheless, a good sensitivity to the mass survives 
and, as we will show, it will be possible to disentangle 
with a good confidence the two cases m v = and 1 
eV. 

To test quantitatively the idea outlined above we 
proceed in two steps: 

i) First we generate a set of synthetic neutrino signals, 
according to some suitable SN model. Neutrinos are 
then propagated from the SN to the detector assuming 
two different SN-Earth distances (10 and 20 kpc) and 
two different mass values m v — and 1 eV. Finally, 
two different energy thresholds (5 and 10 MeV) are 
used for the detection. The result consists of several 
neutrino samples, that hopefully would not differ too 
much from a real SN signal. 

ii) The signals are then analyzed with the main goal of 
disentangling with sufficient confidence the two cases 
m v — and 1 eV. Only the SN-Earth distance is as- 
sumed to be known (this information could be ob- 



tained directly from optical observations or, with an 
uncertainty of the order of 25%, from a comparison 
of the measured total energy with a theoretical esti- 
mation of the binding energy released [13 )• Other 
quantities, like the spectral functions and the details 
of the time evolution of the neutrino flux, are inferred 
directly from the data. 

Generation of the neutrino samples. The time evo- 
lution of the neutrino luminosity and average energy 
can be obtained from the results of SN explosions sim- 
ulations 0, [23, 0, • In the present analysis we 
use the results of Woosley et al. 20]. This model 
is characterized by a rather hard neutrino spectrum 
{(E Pc ) > 20 MeV) that results in a large number of 
detected events (8,800 in SK). Still this is a conser- 
vative choice, since it also implies a depletion in the 
number of low energy neutrinos with a correspond- 
ing loss of informations on the neutrino mass. Inside 
the protoneutron star neutrinos are in thermal equilib- 
rium, and therefore it is reasonable to assume that the 
gross features of their energy spectrum after emission 
can still be described by a Fermi-Dirac distribution: 
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where F 2 (0, if) is defined in (0J and a 'pinching' factor 
77(f) has been introduced to simulate spectral distor- 
tions jl7j. From now on, quantities with a hat (fj, T, 
. . . ) represent inputs to the Monte Carlo (MC), while 
quantities without a hat will refer to results of data 
fitting. We define the Fermi functions and generalized 
Fermi integrals as: 
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Using the time evolution of the average neutrino en- 
ergy E u {i) as given in fig. 3 of ref. |2J| and taking 
for simplicity 77 = 3 0] constant, we compute the ef- 
fective temperature that, together with 77, determines 
the neutrino spectrum at the source: 



T(f) = k 3 E v (t) , 



(•5) 



where k 3 = F 2 (0, 3)/F 3 (0, 3) w 4. Given that the av- 
erage energy is only mildly dependent on time |20| , we 
model the evolution of the neutrino flux $(f) simply 
by taking it prop ortional to the luminosity as given in 
fig. 2 of ref. Hof. 

In the SK detector, z7 e 's are detected through the 
inverse j3 decay reaction v e p — > e + n that has the 
energy threshold -E roa ct = m e + Am„ p ~ 1.8 MeV, 
where m e is the electron mass and A.m np ~ 1.3 MeV 



3 



is the neutron to proton mass difference. The ef- 
fect of the detection cross section a{E v ) 23] is taken 
into account from the beginning by using the function 
S(E„;f{t),r)) x $(t) x a{E v ) in the MC generator. 
Each SN neutrino is labeled by its emission time t u 
and by its energy E v . The corresponding detected 
positron is also identified by an energy/time pair of 
values (E e ,t e ). We generate E e according to a Gaus- 
sian distribution with central value E e = E v — Am np 
and variance a = 0.15 -^10 MeV/i? e that corresponds 
to the SK energy resolution |24j. Since the time res- 
olution of the SK detector is very precise, no error 
is assigned to t e . For massive neutrinos, we rede- 
fine t e = t v + M v by including the appropriate de- 
lay Af„ = L ml/2 El. A fixed number of 8,800 en- 
ergy/time pairs (Ef, f|) is generated in each run. This 
corresponds to the expected number of u e interactions 
within the SK fiducial volume for the spectrum of the 
model in [2(J and a SN at 10 kpc. For L = 20 kpc this 
number is resampled down by a factor of four. Next, 
positrons with energies below the detection threshold 
(Etr = 5 or 10 Mev) are discarded, and as a last step 
the origin of the time axis is set in coincidence with the 
first positron detected (tf = 0) and each f| is accord- 
ingly redefined. For each set of values (E tI ,L,m v ) 
forty different samples are generated in this way. 

Analysis of the neutrino signal. We analyze the 
neutrino signals by means of the likelihood function 



L = S(e; T(t), 77(f)) x $(i + St; b, d, f) x a{e 
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The energy e = E e +Am np of each neutrino is inferred 
from the positron energy. The spectral function S is 
assumed of the form (J2J with a time dependent effec- 
tive temperature T(t) and pinching factor 77(f) fitted 
from the data. &(t + St; b, d, f) describes the time evo- 
lution of the neutrino flux, and the four parameters 
St ,6 ,d and / account for its location on the time axis 
and detailed shape. Finally, cr(e) is the neutrino cross 
section. 

The spectral functions. The two spectral functions 
T(t) and 77(f) are determined by fitting the first 
and second momentum of the energy distribution to 
the mean energy E v and mean squared energy E% 
of the incoming neutrino flux (prior to detection). 
Given a set of n neutrinos with measured energies 
(ei, £2, ■ • ■ -e n ), and assuming for the cross section the 
approximate form o~(e) oc e 2 , we have: 



E „ E, ^ 
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From (JTJ we see that the typical neutrino delays are 
at most of the order of few milliseconds. This is much 
shorter than the time scale over which sizable vari- 
ations of the spectral functions are likely to occur 
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FIG. 1: Comparison of the Monte Carlo spectral temper- 
ature T{i) and pinching parameter i) (thick lines), with 
T(t) and 77(f) fitted for 40 samples (L = 10 kpc, m v = 0). 

Therefore, by binning over sufficiently 
large time windows we can ensure that the determi- 
nation of T and 77 for each bin will not be affected 
by the delays. We use a set of 20 windows of size 
increasing from 50 ms to 2 s, distributed over a sig- 
nal duration 20 s. After this time the neutrino flux is 
assumed to be too low to provide additional informa- 
tions. To reduce further the possible effects of mass 
induced time delays, only neutrinos with energy larger 
than E* = 15MeV are used. For each time window 
[*ai t a +i\ we perform a least square fit to the following 
quantities: 
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where the functions F n (y, 77) defined in flf. depend on 
the temperature through y = E* /T. This yields the 
'best' values T(t ) and ?7(f a ) at i a = (t a +t a+1 )/2. Fi- 
nally, in order to obtain two smooth continuous func- 
tions, the points T(t a ) and rj(i a ) are interpolated with 
two half integer power polynomials Pt.^ ~ J2i c '^ 2 
with I = 0,1, 2,... 10 for T and / = 0,1,3,5 for 77. 
In fig. ^ the functions T(t) and 77 = 3 used in the 
MC (thick lines) are compared with forty fits to neu- 
trino samples generated with m v — and L = 10 kpc. 
These results remain unchanged for 771^ = 1 eV, while 
for L = 20 kpc, due to the reduced statistics, the fits 
show a somewhat wider dispersion. 

The neutrino flux. Even if the details of the neu- 
trino flux evolution with time are not known, its gross 
features can be predicted on rather solid theoretical 
grounds. It is expected that a sharp exponential rise, 
with a time scale of tenths of milliseconds, is followed 
by a power law decay, with a time scale of several 
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seconds. We model this behavior by means of a para- 
metric analytical function in which, roughly speaking, 
two parameters describe the rising and decaying rates, 
and a third one accounts for the transition point: 

-jfr^^i^" lZ%. ^ 

This function can reproduce reasonably well the re- 
sults of different models [H El The two 
exponents m and n are fixed to suitable integer values 
by means of a preliminary rough fitting procedure to 
the data, and then are held constant throughout the 
analysis of all the samples. For the model in [i(J we 
use n — 8 and m = 2. Finally, since the origin of times 
was arbitrarily set in coincidence with the detection of 
the first neutrino while <!> vanishes at t = 0, a fourth 
parameter St is required to let <&(£ + St) freely shift 
along the time axis. 

The likelihood analysis. The minimization of the 
negative log-likelihood is carried out by means of the 
package MINUIT 25]. To avoid double minimums, we 
minimize with respect to the square of the neutrino 
mass. Given a value of the time delay of each 
neutrino is computed according to its energy ej, and 
subtracted from its arrival time The log-likelihood 
for the new array of times is then evaluated, and min- 
imized with respect to the other parameters b ,d, f 
and St. This proceeds until the absolute minimum 
is found in the full five dimensional parameter space. 
There is a subtlety related to the cases when, espe- 
cially for large test masses, a neutrino is migrated to 
an early time where $ = 0. The problem is not just a 
numerical one of logarithm overflow. Due to the un- 
certainty in the energy measurement, the firsts neu- 
trinos detected can end up in such a position without 
necessarily implying that the corresponding distribu- 
tion has vanishing probability. To account for this, for 
the relevant neutrinos the energy uncertainty is con- 
verted into an uncertainty in the new time position, 
and the corresponding contribution to the likelihood 
is evaluated by convolving C{t) in (fTUf) with a Gaus- 
sian of the appropriated width. In fig. [21 the results 
of the fitted fluxes for 40 samples with L — lOkpc, 
E tI = 5 MeV and m„ = eV are compared with the 
function used in the MC generator (thick line). 

Results. Our results are summarized in fig. [3] and in 
table 1. Fig. |31 depicts the results of the mass fits for 
two sets of 40 samples with E tr = 5 MeV, L = 10 kpc, 
rh v = (circles) and rh v = 1 eV (squares) . The data 
points have been ordered with increasing value of the 
best fit mass to put in evidence the difference be- 
tween the two cases. Error bars correspond to 95% 
c.L, where the lower (upper) limit is computed by in- 
tegrating the likelihood from = oo (— oo) until 
reaching the 95% of the area, while minimizing with 



*(t) 




0.1 sec. 0.2 0.3 

FIG. 2: Comparison between the Monte Carlo flux func- 
tion (thick line), and the result of the likelihood anal- 
ysis for 40 samples (L — 10 kpc, Etr = 5 MeV, m„ = 0). 

respect to the other parameters. We have checked 
that the limits do not change much if the integration 
is restricted to the physical region rr? v > 0. 

Table 1 summarizes some results of the mass fit- 
ting. For each set of parameters (-Etr) L, m„) we have 
analyzed 40 samples. The first raw gives the percent- 
age of times in which the 95% c.L lower limit m l v is 
larger than the input mass rh v . The second raw refers 
to the cases when the upper limit m" is smaller than 
rh v . Number in parentheses correspond to L = 20 kpc. 
These figures characterize the percentage of 'failures' 
of the method, that therefore appears to be reliable in 
about 90%-95% of the cases. The third row gives the 
percentage of times when = is excluded at 95% 
c.L when the signal is generated with m v — 1 eV. This 
characterizes the power of the method for excluding a 
massless neutrino. We see that in the most favorable 
case (Etr = 5 MeV, L = 10 kpc) the method is suc- 
cessful in more than 50% of the cases. The following 
three raws give the average over the 40 samples of the 
mass square best fit and of the 95% c.L lower and up- 
per limits (only for L = 10 kpc). These last figures are 
just intended to give an idea of the quality of the fits 
(they would be fully meaningful only if 40 SN could 
be observed). 

From the results in the table it is apparent that 
low energy neutrinos are crucial for the sensitivity of 
the method, and therefore a low detection threshold is 
very important. With E tI = 5 MeV the massless case 
is excluded in about 50% of the cases, while this drops 
to 25% when E tI = 10 MeV. Also, with the higher en- 
ergy treshold the fluctuations of the results over the 
40 runs is doubled (last three rows). Unfortunately, 
there could be a dangerous background in the energy 
range between 5 and 10 MeV, represented by pho- 
tons originating from neutral current reactions off 16 
mainly produced by the more energetic fi and r neu- 
trinos |26(. Softer neutrino spectra, like for example 
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FIG. 3: Fitted values of the mass and 95% c.l. error bars 
for two sets of 40 analysis with L = lOkpc and E tr = 
5MeV, for m„ = (circles) and ra„ = leV (squares). 
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0.6 ±0.5 
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TABLE I: Results of the analysis for the two neutrino 
masses m„ = and 1 eV, two energy thresholds of 5 and 
lOMeV and for two SN distances of lOkpc and, in paren- 
thesis, 20kpc. The first (second) row give the percentage 
of times in which the 95% c.l. lower (upper) limit m v (m™) 
is larger (smaller) than rh v . The third row lists how many 
times m„ — can be excluded (95% c.l.) when m„ = 1 eV. 
The following three rows give, for L = lOkpc, the mass 
square best fit and the 95% lower and upper limits aver- 
aged over 40 samples. 

the ones predicted by the model of Totani et al. |2l| , 
would have the double benefit of drastically reducing 
this background, while at the same time raising the 
number of v e emitted in the 5-10 MeV energy range 
from about 7% of the present analysis to about 20%. 
In this case a better sensitivity to the mass could be 
expected. We should also mention that very recently 
it has been suggested that water Cerencov detectors 
could be modified to allow tagging of the inverse f3 
decay neutrons This would eliminate the back- 

ground from neutral current reactions and allow for 
lower thresholds. 

Numerical spectrum and neutrino oscillations. The 
previous results have been derived relying on two 
main simplifying approximations for the neutrinos en- 
ergy spectrum: 1 ) the neutrino energies have been 
generated assuming the 'pinched' Fermi-Dirac spec- 



trum given in eq. @ and fitted, as was described 
above, with a similar two-parameters energy distri- 
bution; 2) no effects of the neutrino oscillations were 
included in the analysis. In order to evaluate to what 
extent the sensitivity of the method could be affected 
by these approximations, we have run a set of simu- 
lations in which the neutrino energies were generated 
according to the shapes of the numerical spectra given 
by Janka and Hillebrandt in ^t|- A time dependent 
energy rescaling of the spectral shapes was introduced 
to reproduce properl y th e time evolution of the mean 
energies as given in |2fJ. We stress that since a two 
parameters Fermi-Dirac distribution fits rather accu- 
rately the spectra obtained from the numerical simula- 
tions, dropping our first simplification does not affect 
sensibly the numerical results. 

Neutrino oscillations can produce a composite spec- 
trum corresponding to an admixture of the origi- 
nal v e spectrum with a harder component due to v x 
(x = h,t) [iU. Clearly, the resulting spectral dis- 
tortions will depend in first place on the size of the 
v e — v x spectral differences. While it is often stated 
that these differences could be quite sizable, and could 
yield up to a factor of two hierarchy between the v x 
and v e average energies 0, |2(1 0, H^, recent and 
more complete analysis of SN neutrino spectra forma- 
tion 29] indicate that this is not the case: the inclu- 
sion of important interaction rates that were neglected 
in previous works yields spectral differences that are 
only of the order of 10% [29j. At this level, identify- 
ing the two components of a mixed spectrum would 
be a difficult task, and could represent a real chal- 
lenge for the study of SN neutrino oscillations in the 
v e — Vx channel. However, for what concerns our anal- 
ysis, this ensures that the results discussed above are 
not affected much by neutrino oscillations. To keep 
on the safe side, we have tested the sensitivity of our 
method to oscillation effects by running a set of sim- 
ulations using the results of Woosley et al. [2(j for 
which, as a consequence of negelcting important neu- 
trino reactions |29j, the spectral differences are ex- 
treme {Ejj^, w 2Eu,.). As we will show, even in this 
(probably unrealistic) case, we find that the loss in 
sensitivity to the neutrino mass is small. We can con- 
clude that the effects of neutrino oscillations do not 
endanger the applicability of our method and neither 
its overall sensitivity to the neutrino mass. 

As is discussed in [2^, depending on the type of 
the neutrino mass hierarchy (normal or inverted) and 
on the size of sin 2 #13 (large > 10~ 4 or very small 
< 10~ 6 ), neutrino oscillations could i) harden the v e 
spectrum through a complete spectral swap with 
(inverted hierarchy, large #13) ii) mix the v e spec- 
trum with a fraction of about sin 2 6*12 ~ 1/4 of 
harder neutrinos (in the other cases). (The region 
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FIG. 4: Comparison between the i? e (upper thick line) and 
j?fj, (lower thick line) fluxes used in the Monte Carlo, the 
resulting v e flux at the detection point (middle thick line), 
and the maximum likelihood fits for 40 neutrino samples 
with a mixed composite spectrum. 
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FIG. 6: Fitted values of the mass and 95% c.l. error bars 
for a set of 40 neutrino samples with mixed composite 
spectra (L = lOkpc, E tr = 5MeV, m„ = (circles), rh v — 
1.2 eV (squares)). 



intermediate amount of mixing.) The first case can 
be studied without modifications in our procedure. 
Clearly, a different spectrum would imply somewhat 
different numerical results; however, this is analogous 
to the unavoidable uncertainty related to the choice 
of the particular SN model since, for example, the v e 
mean energies that have been used in the present anal- 
ysis 113] are quite close to the mean energies of the 
model of Totani et al. 21] . The second case is more 
interesting since, for large spectral differences, fitting 
a composite spectrum with just one effective spectral 
temperature and one pinching parameter could de- 
grade somewhat the sensitivity. 

We have carried out an analysis of forty neutrino 
samples with a mixed composite spectrum as would 
result from a normal mass hierarchy, sin 2 d\i w 1/4 
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FIG. 5: Comparison between the v R (lower thick line) and 
Vfi (upper thick line) mean energies of the numerical spec- 
tra used in the Monte Carlo, the mean energy of the mixed 
composite spectrum (middle thick line), and the fitted val- 
ues of _Ej7 c (t) for 40 neutrino samples. 



and large $13. The original v e and D x fluxes taken 
from |2C| are depicted in fig. 0] and compared with 
the composite v e flux at the detection point, as well 
as with our maximum likelihood fits. The relative 
normalization of the fluxes was computed assuming 
flavor equipartition of the integrated luminosities and 
the time dependent average energies given in j^. A 
similar comparison between the mean energies of the 
two spectral components, the average energy of the 
mixed spectrum and the average energies obtained 
from Fermi-Dirac distributions with fitted spectral pa- 
rameters T(t) and fj(t) is depicted in fig. El 

For the mass fits we have taken E tr = 5 MeV and 
L = 10 kpc. Since we are essentially interested in the 
loss of sensitivity to the neutrino mass with respect to 
the non-oscillation case, we have been searching for a 
mass value that can reproduce results similar to those 
of the band plot in fig. |3| (namely, for a signal gen- 
erated with rh v ^ we require the massless case to 
be excluded at 95% c.l. in at least 50% of the runs). 
As is shown in fig. H3 our requirement is fullfilled for 
rh u = 1.2 eV, to be compared with m„ = 1.0 eV for 
the non-oscillation case. We conclude that even in the 
unrealistic case of extremely large spectral differences 
between the components of a mixed spectrum, fitting 
the data with a single two parameters Fermi-Dirac dis- 
tribution does not degrade much the numerical results 
for the mass estimate. Given that the most recent re- 
sults suggest that the v e and spectra are in fact 
not very different [2^, the use of more complicated 
bimodal energy distributions to improve the mass fits 
is probably not justified. 

Before concluding, a few remarks are in order. In 
this work we have not carried out any deep study 
aimed to optimize the single fits to each different neu- 
trino sample (windows size, interpolating functions 
Pt,ti , specific flux function parameters) . A large num- 
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ber of different sets of data have been analized through 
the very same procedure, in order to collect enough 
informations about the sensitivity of the method in 
a reasonable amount of time. It is clear to us that 
optimizing the overall procedure in order to analze a 
specific sample (as would certainly be the case with a 
signal from a real SN) can improve the sensitivity and 
shrink somewhat the uncertainties on . 

Besides the effects of neutrino oscillations that were 
briefly analyzed in the last paragraph, a few other is- 
sues could deserve further investigation. For example, 
assuming that the SK signal could be combined with 
negligible uncertainty on the absolute timing with the 
signal detected at KamLAND, the much better en- 
ergy resolution and the lower threshold of this last 
detector could enhance the sensitivity to the neutrino 
mass. It also remains to see what sensitivity could 
be achieved with a statistics one order of magnitude 
larger, as would be available with the megaton neu- 
trino detectors presently under study 
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